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The iterated Crank-Nicholson method has become a pop- 
ular algorithm in numerical relativity. We show that one 
should carry out exactly two iterations and no more. While 
the limit of an infinite number of iterations is the standard 
Crank-Nicholson method, it can in fact be worse to do more 
than two iterations, and it never helps. We explain how this 
paradoxical result arises. 
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I. INTRODUCTION 

There is currently a large worldwide effort underway 
attempting to solve Einstein's equations numerically for 
astrophysically interesting scenarios. The problem is ex- 
tremely challenging technically. Among the difficulties 
is that of finding a finite-difference scheme that allows a 
stable time evolution of the system. It is well-known that 
implicit differencing schemes tend to be stable. However, 
the difficulty of solving the resulting implicit algebraic 
equations, especially in three spatial dimensions, has led 
most researchers to stay with explicit methods and their 
potential instabilities. 

Several years ago, Choptuik proposed solving the im- 
plicit Crank-Nicholson scheme by iteration. This would 
effectively turn it into an explicit scheme, but hopefully 
by iterating until some convergence criterion was met one 
would preserve the good stability properties of Crank- 
Nicholson. The iterated Crank-Nicholson scheme has 
subsequently become one of the standard methods used 
in numerical relativity. 

In this note, we point out that when using iterated 
Crank-Nicholson, one should do exactly two iterations 
and no more. While the limit of an infinite number of 
iterations is the implicit Crank-Nicholson method, it can 
in fact be worse to do more than two iterations, and it 
never helps. 



II. ITERATED CRANK-NICHOLSON 



are similar.) A simple first-order accurate differencing 
scheme is FTCS (Forward Time Centered Space): 
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Here n labels the time levels and j the spatial grid points. 

It is a standard textbook result that this scheme is 
unconditionally unstable. One sees this with a von Neu- 
mann stability analysis: Put 
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and find that the amplification factor ^ is 

^ — 1 + ia sin kAx, 

where a = At/ Ax. Since > I for any choice of 
the method is unconditionally unstable. 

Backwards differencing gives a stable scheme: 



(3) 



(4) 



a. 



.n+l 



.,n+l 



,n+l 



for which 



At 



2Aa 



1 



ikAd 



(5) 



(6) 



Now < 1 for any choice of a, and so the method is 
unconditionally stable. 

The Crank-Nicholson scheme is a second-order accu- 
rate method obtained by averaging equations and (||) . 
Now one finds 
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Since = 1, the method is stable. It is the presence 
of the quantities on the right hand side of equation 
(^ that makes the method implicit. 

The first iteration of iterated Crank-Nicholson starts 
by calculating an intermediate variable ^^^tt using equa- 
tion (0): 



To understand this paradoxical result, consider differ- 
encing the simple advective equation 



du du 
dt dx 
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(Many equations in numerical relativity are general- 
izations of this form, and the differencing techniques 
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Then another intermediate variable '•^•'u is formed by av- 
eraging: 



l((l)^n+l+^n). 



(9) 
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Finally the timestep is completed by using equation 
again with u on the right-hand side: 
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(Iterated Crank-Nicholson can alternatively be imple- 
mented by averaging the right-hand side of equation (|^). 
For linear equations, this is completely equivalent.) 

Iterated Crank-Nicholson with two iterations is carried 
out in the same way. After steps (||) and (0), we calculate 



(2)~n+l 



(1) «+l/2 _ (l)-«+l/2 



At 2Aa; 



(11) 
(12) 



Then the final step is computed analogously to equation 
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Any number of iterations can be carried out in the same 
way. 

Now consider the stability of these iterated schemes. 
If we define (3 = (a/2) sin/cAa;, and caU the FTCS 
scheme (||) the zeroth-order method, then direct calcula- 
tion shows that the amplification factors are 



(•'^e = 1 + 2i/3, 

(1) ^ = l + 2^/3-2/3^ 

(2) ^ = 1 + 2^/3-2/3^-2^/33, 

(3) ^ = 1 + 2i/3 - 2/32 - 2i(3^ + 2/3^ 



(14) 
(15) 
(16) 
(17) 



and so on. As one would expect, these are exactly the 
same values one gets by expanding equation (|^) in powers 
of (3 and truncating at the appropriate point. 

To check stability, compute for each of these ex- 
pressions. You find an alternating pattern. Levels and 
1 are unstable; levels 2 and 3 are stable provided < 1; 
levels 4 and 5 are unstable; levels 6 and 7 are stable 
provided /3^ < 1; and so on. Since the stability require- 
ment must hold for all wave numbers fc, it translates into 
a^/4 < 1, or At < 2Aa;. This is just the Courant condi- 
tion (the 2 occurs because of the 2 in eqn. |^). 

Now we see the resolution of the paradox: While the 
magnitude of the amplification factor for iterated Crank- 
Nicholson does approach 1 as the number of iterations 
becomes infinite, the convergence is not monotonic. The 
magnitude oscillates above and below 1 with ever de- 
creasing oscillations. All the cases above 1 are unstable, 
although the instability might be very slowly growing for 
a large number of iterations. 

The accuracy of the scheme is determined by the trun- 
cation error. This remains second order in At and Aa; 
from the first iteration on. Doing more iterations changes 
the stability behavior, but not the accuracy. Since the 



smallest number of iterations for which the method is 
stable is two, there is no point in carrying out more iter- 
ations than this. 

Note that there was nothing special about using the 
advective equation (^) for this analysis. Similar behavior 
is found for the wave equation, written in first-order form 
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with the standard centered difference formula for the sec- 
ond derivative term. One recovers the usual Courant 
condition (without the factor of 2) for the stable cases. 
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